Regular modes in rotating stars 
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Despite more and more observational data, stellar acoustic oscillation modes are not well under- 
stood as soon as rotation cannot be treated perturbatively. In a way similar to semiclassical theory 
in quantum physics, we use acoustic ray dynamics to build an asymptotic theory for the subset of 
regular modes which are the easiest to observe and identify. Comparisons with 2D numerical simula- 
tions of oscillations in polytropic stars show that both the frequency and amplitude distributions of 
these modes can accurately be described by an asymptotic theory for almost all rotation rates. The 
spectra are mainly characterized by two quantum numbers; their extraction from observed spectra 
should enable one to obtain information about stellar interiors. 
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Stars being far-away objects, the types of information 
that can be obtained from them are necessarily limited. 
One of the most important corresponds to luminosity 
variations, which can reflect the passing of a planet or 
intrinsic modulations in the light emitted by the star. In 
particular, the domain of asteroseismology studies stellar 
oscillation modes, which create periodic variations of the 
luminosity which can be detected [l| . For the Sun, these 
modes have been theoretically constructed and success- 
fully compared with observations, leading to detailed in- 
formation on the Sun's internal structure. However, this 
theory requires the star to be nearly spherically symmet- 
ric, an assumption clearly violated for rapidly rotating 
stars 0]. With the launch of the recent space missions 
COROT and Kepler [l], oscillation spectra of rapidly 
rotating stars are observed with great accuracy. This 
concerns mainly the stars more massive than the Sun 
that belong to the main sequence of the Hertzsprung- 
Russel diagram. In order to access their internal struc- 
ture through a seismic diagnostic, it is thus essential 
to understand the oscillation spectra of rapidly rotating 
stars 

Accurate computations of acoustic modes fully taking 
into account the effects of rotation on stellar oscillations 
have only recently been performed for rotating stars [3] 
(an example is shown in Fig. [T]). Such stationary pat- 
terns of acoustic waves can be described asymptotically 
through their short- wavelength limit, in the same way 
as classical trajectories can describe quantum or electro- 
magnetic waves in this limit Q . These acoustic rays obey 
Hamiltonian equations of motion. In [5], their dynamics 
was investigated for a polytropic stellar model, showing 
that the tools from the fields of classical and quantum 
chaos enable us to understand the behavior of modes in 
rapidly rotating stars. Indeed, for increasing rotation 
rates, the dynamics undergoes a transition from an inte- 




FIG. 1: (Color online) Pressure amplitude Py/d/po on a 
meridian plane for a polytropic model of stars, with d the 
distance to the rotation axis and po the equilibrium density. 
The mode shown corresponds to n — 46, 1 = 1 and m = 
at a rotation rate of = 0.783, with Q, K = (GM/Rl q ) 1/2 
being the limiting rotation rate for which the centrifugal ac- 
celeration equals the gravity at the equator, M the stellar 
mass and R eq the equatorial radius. Colors/grayness denote 
pressure amplitude, from red/gray (maximum positive value) 
to blue/black (minimum negative value) through white (zero 
value) . The thick black line on the right is the central periodic 
orbit 7 of the island. 



grable to a mixed system, where chaotic and stable zones 
coexist in phase space. The asymptotic theory built for 
slowly rotating stars, which does not take these effects 
into account, cannot thus be applied at high rotation 
rates. In the latter regime, it was shown that the spec- 
trum of acoustic oscillations can be divided into several 
subspectra corresponding to regular and chaotic zones in 
phase space in a way similar to what happens in quantum 
chaos systems [f|. 
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FIG. 2: (Color online) Surfaces of section at rotations = 
0.224 (top) and j£- = 0.589 (bottom) for m = 0. Each dot 
represents the crossing of an acoustic ray with the equatorial 
half-plane, r being the radial coordinate and k r the associ- 
ated momentum. Orange/light gray denote a chaotic ray, 
green/dark gray a whispering gallery ray, blue/black a stable 
island ray (see text). Upper insets are close-ups of the main 
stable island. Lower inset in the top figure shows a close-up 
of the main island at = 0.262, just after the bifurcation. 



As already demonstrated for the Sun and solar-type 
stars, a quantitative asymptotic theory is crucial to ex- 
tract information from the observed spectra as it links 
the behavior of oscillation modes to physical properties of 
the star Q ■ Furthermore, the complexity of the observed 
spectra usually requires prior knowledge of an asymp- 
totic theory in order to correctly identify the frequency 
peaks in the data with specific oscillation modes. In this 
paper, we present for the first time such a quantitative 
theory at almost all rotation rates for a specific subset 
of modes, which should be among the easiest to obtain 
from observations. Indeed, we focus on a series of modes 
centered around the largest stable island of the system, 
and systematically build them using the parabolic equa- 
tion method 0. This method was successfully applied 
to light in dielectric cavities @, electronic resonators in 
a magnetic field [l(| and quantum chaos systems fTlj j . 
The results of the asymptotic modes are then compared 
with numerical computation of oscillations in a poly- 
tropic star, showing that the theory correctly describes 
the modes even in the bounded frequency range of stellar 
oscillations. 



The study of ray dynamics in [B| showed that for a wide 
range of rotation rates, three main types of phase space 
zones with different dynamics can be defined (see FigJ2J : 
1) regular structures built around stable periodic orbits 
(stable islands); 2) whispering gallery rays close to the 
surface; 3) chaotic zones with ergodic rays in the remain- 
ing parts of phase space. Figure [2] shows that for m = 
(axisymmetric modes) the main island undergoes a bifur- 
cation from one island centered on the rotation axis to 
two islands which move away from the rotation axis as the 
rotation rate increases. Each phase space region gives rise 
to a well-defined subspectrum of modes localized inside 
the region. The whispering gallery modes are essentially 
unobservable in real stars since the disk-average leads to 
a very small contribution in observed spectra |5| . Chaotic 
modes can have visible contributions, but the associated 
spectra, although they can display well-defined statisti- 
cal properties, cannot be described by a few quantum 
numbers. In contrast, the stable island modes give rise 
to very regular sequences of frequencies described by a 
few parameters which can be potentially extracted from 
observed spectra. These modes (an example is shown on 
Fig. 1) can be characterized on a meridian plane by the 
number n of nodes along the central periodic orbit 7 and 
the number t of nodes transverse to 7. 

To describe asymptotically these island modes, we 
start from the equation of acoustic waves in stars. We 
neglect the Coriolis force, which is known to be negligi- 
ble in the high-frequency regime since the Coriolis force 
time scale (l/(2f2)) is much longer than the mode period 
[H . We also neglect the perturbations of the gravitational 
potential, since they are produced by the density fluctu- 
ations and tend to cancel out as the number of nodes 
of the density distribution increases for high-frequency 
modes, as has been numerically checked for non-rotating 
stars 12]. Finally, we use the adiabatic approximation 
which is known to be a very good approximation to com- 
pute frequency modes in non-rotating stars. Indeed, it 
is accurate enough to interpret the solar acoustic modes 
despite the fact that these frequencies are determined to 
high accuracy [;7] . In the linear approximation, this gives 
rise to a Hclmholtz-type equation. Using the cylindri- 
cal symmetry of the system with respect to the rotation 
axis, one can rewrite this equation as a two-dimensional 
problem 



c 2 (m 2 



d 2 



(1) 



Here $ m is the mode amplitude scaled by the square root 
of the distance d to the rotation axis, c s is the sound ve- 
locity (which depends on the location inside the star), to 
is the frequency of the mode, and u> c is the cut-off fre- 
quency whose sharp increase in the outermost layers of 
the star leads to the reflection of acoustic waves and thus 
the formation of modes in a bounded frequency range 
w < (^c)max- The integer m is the quantum number 
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corresponding to the quantization of angular momentum 
along the rotation axis. To construct the island modes 
centered on the stable periodic orbit 7 of length L 7 , we 
rewrite Eq. ^ in the coordinates (s, £) centered on 7, 
with s the coordinate along 7 and £ the transverse co- 
ordinate. The parabolic equation method [1] assumes 
that the solutions have a longitudinal scaling in 1/u/ and 
a transverse scaling in 1/y/Uj. We introduce the WKB 
ansatz: 



$ m (s,£) = exp(iuT)U m (s,£,u}), 



(2) 



with dr = ds/c s , using the renormalized sound velocity 



defined by c s = <r 8 ur / (w 
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An expansion in powers of to yields a series of equa- 
tions; keeping terms of order uj and introducing the vari- 
able v = y/oJS, and the function V n , 
the parabolic equation: 



U m /y/c s (s) give 



d 2 V m 
dv" 1 



2i 



1 dV m 

c s (s) ds 



K{s)v 2 V m 



(3) 



where K(s) = j^^i 



The equation in the trans- 
verse coordinate v is similar to the harmonic oscillator 
in quantum mechanics, with an additional term depend- 
ing on the longitudinal coordinate. The ground state is 



of the form: = A(s) exp 



equations: 



dT(s) 



•r(s) 2 



and obeys the two 



3 (s) ds 



■T(s) 2 + K{s) = Oand 



dA(s) _ 



A(s 



dz(s) 



_M£) r ( S ). Defining z(s) through T(s) = da 

implies that A(s) = 1/ 'y/z(s) and z(s) should satisfy the 
following system of equations: 



1 dz 
c s (s) ds 



P , - 



1 dp 

c s (s) ds 



-K(s)z 



(4) 



Equations (UJ) are periodic in the variable s with period 
Ly, and according to Floquet theory there exists an op- 
erator describing the evolution over one period. To con- 
struct it, one uses the fact that system (Q| corresponds to 
the Hamilton equations associated with the Hamiltonian 



H= £r 



K(s 



-sj- 1 — 2 — ■ ^ can ^ e shown that the same equa- 
tions (with z and p real) describe the acoustic ray in the 
vicinity of the central periodic orbit (via a normal form 
approximation), z being the transverse deviation from 7 
and p the associated momentum. We linearize the motion 
around the periodic orbit and construct the monodromy 
matrix which describes this linearized motion from one 
point to its image after one period: 



Z(S + L-y) 



M 



z(s) 
p(s) 



(5) 



For the mode to be univalued, V„ t should be the same 
after one period up to a global phase and thus z and 



p should correspond to an eigenvector of M. As 7 is 
stable, the matrix M is conjugate to a rotation matrix 
and has two eigenvalues e ±m , with a in [0, 2n[. The 
corresponding eigenvectors are complex conjugate, only 
one of them giving the physical solution exponentially 
decreasing at large v. 

The modes of higher frequency can be constructed 
as for the harmonic oscillator from using standard 
methods from quantum mechanics; the result, up to 
a normalization constant, is equivalent to multiplying 
V^{s,v) = z~ 2 exp [i^ 2 ] by a function containing the 
Hermite polynomials of order I noted Hf. 



exp 



(6) 



Again, for the mode to be univalued, the global 
phase accumulated after one period should be a multi- 
ple of 2-7T. This phase is exp exp (—i(2wN r + at)/ 2) 

exp(-i(27riV r + a)£) exp (iuj §^ The first two 

phases correspond to the so-called Maslov indices [1, [l3| 
and count the number of caustics encountered in the lon- 
gitudinal and transverse motions. The number N r keeps 
track of the number of times the trajectory solution of 
Eq. |4]) makes a complete rotation around 7 in phase 
space, and can be evaluated from ray simulations. This 
implies: 



If c„ 



27r(rc+-) 



'A 



(2irN r + a) 



(7) 

The regular subspectrum is thus essentially described 
by two quantities, Sn = -pjr and 81 = 2 ^ V dt"° (which 

depend on m). This corresponds to the empirical for- 
mula found in 3} for m = from numerical simulations. 
The quantities Sn and 51 probe the sound velocity along 
the path of the periodic orbit and its transverse deriva- 
tives. Indeed, an explicit expression of a in terms of such 
transverse derivatives can be derived Equation ([7]) 
is valid asymptotically for n large and i <?C n. As observ- 
able modes in real stars cannot be too high in frequency, 
we have checked numerically the validity of this formula 
for moderately high values of n. 

In Fig|3]we plot the numerically computed Sn and St 
values vs. the theoretical ones for a large range of ro- 
tation rates. We restrict ourselves to the case \m\ < 1 
which is the most common in observational data. Nu- 
merical modes were obtained using a code that computes 
adiabatic modes of rotating polytropic stars as in Q , and 
selecting the island modes through their phase space lo- 
cations. Theoretical values were obtained from the the- 
ory explained above, estimating the monodromy matrix 
entries by following classical trajectories in the vicinity 
of the periodic orbit, using the fact previously noted that 
Eq. (0| describes the deviation of a nearby classical tra- 
jectory from the central orbit 7. We checked that the 
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FIG. 3: (Color online) Comparison between actual regulari- 
ties of regular modes and theoretical predictions for m — 
and different values of SI/Qk (ui p = y/GM/Rp with R p the 
polar radius). Circles: 5n, triangles: 51, orange/light gray: 
theory, blue/dark gray: numerical results. Numerical results 
correspond to different sets of values of Sn and 51, with n 
between 42 and 51 and l = Q,l. As m = 0, only modes sym- 
metric with respect to the rotation axis should be retained, 
and thus the theoretical value of 51 is multiplied by two be- 
low the bifurcation point. Close to fl = 0, the separation 
of numerical values in two groups corresponds to the small 
separation in Tassoul's theory p|. Inset: Same for m = 1. 



results were not sensitive to the choice of the trajectory 
inside the island. The results of Figf3] show that a good 
agreement exists between numerical and theoretical reg- 
ularities, except close to ft — where Tassoul's asymp- 
totic theory applies @- For Sn, the agreement is good 
over the whole range of rotation. For 6£, the agreement 
is good at large and low rotation, but degrades in the 
range [0.25, 0.35] for m = 0. We attribute this discrep- 
ancy to the fact that, as seen in Figf2J the periodic orbit 
of the main stability island undergoes a bifurcation in this 
range, from one stable central orbit to two stable orbits 
on each side and a central unstable one. It is known that 
in such a case, the normal form approximation for the 
classical motion which is used in the parabolic equation 
method should be modified by different uniform formulas 
[3| ■ Thus in the vicinity of the bifurcation the method 
is expected not to give accurate results. This picture is 
confirmed by the inset of Figj3l which shows that in the 
case m — 1, where there is no such bifurcation, agree- 
ment is good for 51 over the whole range of f2 values. 
We note that other bifurcations are present in the sys- 
tem which create additional stable and unstable orbits in 
the vicinity of the central one, but they do not seem to 
affect the results for the relatively low-frequency modes 
we consider. We note also that Eq. (JTJ) predicts degen- 
eracies at rational values of a/V. These degeneracies can 
be avoided crossings or true degeneracies if the modes 
belong to different symmetry classes. We have checked 
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FIG. 4: (Color online) Amplitude distributions (real part of 
3> m ) on the equator for a theoretical and a numerical mode. 
Blue dashed line: numerical mode (same as in Fig[l}; orange 
continuous line: theoretical mode, from Eq. 



that it actually enables us to predict such occurrences. 
We note that while the theory neglects the Coriolis force 
and perturbations of the gravitational potential, the nu- 
merical modes were computed taking into account both 
effects. The good agreement seen in Fig J3] confirms that 
these processes can safely be neglected in this regime. We 
also remark that recent analysis of numerically computed 
modes in realistic, non-poly tropic, differentially rotating 
stellar models show the emergence of formulas similar to 
Eq. (|7J| for specific subsets of modes [l5| • 

Not only does the parabolic equation method give the 
frequencies of the modes, but it also yields their ampli- 
tude distribution. Indeed, the eigenvector of the mon- 
odromy matrix gives T(s), which enables us to construct 
an approximation of the mode itself using Eq. ([6]). Com- 
parisons between theoretical and numerical modes show 
that the modes are well approximated by the theory (see 
an example in Fig|4]), although sometimes small oscilla- 
tions due to interference between different modes are not 
well reproduced by the theory. 

In conclusion, we have shown that the parabolic equa- 
tion method enables us to build an asymptotic theory 
for the most visible of the regular acoustic modes of a 
star rotating at arbitrary rotation rates except for very 
slow rotation, where Tassoul's theory Q already applies. 
Comparisons with numerical computations of oscillations 
in a stellar model show that the asymptotic theory gives 
a good description of the frequency differences and am- 
plitude distributions, except for to = at a specific ro- 
tation rate where a bifurcation takes place and a more 
refined theory is needed. The spacings Sn and 5t which 
describe the frequency distribution of this type of modes 
can be expressed in terms of internal characteristics of 
the star. Our results should enable one to use data from 
recent space missions such as COROT and Kepler to ex- 
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tract information about the observed stars and use this 
information to build more accurate stellar models. 
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